Oscillatory Motion of Rising Bubbles in Wormlike Micellar Fluids 
with Different Microstructures 
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Previous observations of the nontransient oscillations of rising bubbles and falling spheres in 
wormlike micellar fluids were limited to a single surfactant system. We present an extensive survey 
of rising bubbles in another system, an aqueous solution of cetylpyridinium chloride and sodium 
salicylate, with and without NaCl, across a range of concentrations and temperatures. Two differ- 
ent types of oscillation are seen in different concentration ranges, each with its own temperature 
dependence. Rheological data allows for the identification of these different hydrodynamic states 
with different fluid microstructures. 

PACS numbers: 47.50.+d, 83.50.Jf, 83.60.Wc 
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Fluids are often broadly categorized as either Newto- 
nian or non-Newtonian, according to whether the Navier- 
Stokes equation does or does not accurately describe the 
fluid's motion. Newtonian fluids consist of molecules 
small enough to be approximated by point masses (such 
as air and water) , while non-Newtonian consist of larger 
structures (such as polymers) . This increased size affords 
greater degrees of freedom which leads to macroscopic 
viscoelasticity 1]. Fluids consisting of self-assembling 
cylindrical (or wormlike) micelles resemble polymeric flu- 
ids on the microscopic scale, with the added feature that 
their length distribution is determined by aggregation ki- 
netics; micelles continually break and reform 0, [jj Q . 
The molecular level physics is clearly more complicated 
for such fluids, and one can reasonably expect this will 
introduce a new set of flow properties to the class of non- 
Newtonian fluids. 

Several novel results have been reported in both vis- 
cometric flows [a, p, B , B , l al and hydrodynamic (non- 
viscometric) flows [l(J, UJi °f various wormlike micel- 
lar fluids. Among these, a falling pendant drop at the 
bottom of a thin filament has been observed to slow to 
a complete stop before the filament suddenly ruptures 
|lfj|; this rupture has been confirmed in other configu- 
rations Also, rising air bubbles and falling solid 
spheres have been observed to oscillate without reaching 
a terminal velocity 0, ^2 • While in polymer solutions 
a rising bubble displays a sharp non-axisymmetric cusp, 
which remains unchanged during a steady rise P, Il4| , in 
wormlike micellar fluids the cusp periodically extends to 
a sharp point, then retracts to a blunt edge (see Fig.l). It 
is likely that these new dynamics are a manifestation of 
the reversible scission reactions of the micelles, though 
the precise mechanism may be quite complicated. An- 
other type of fluid is known to include scission-like reac- 
tions at its microscale - associating polymer solutions - 
and falling sphere oscillations are reported in that system 
as well [lj [l6|. 

In this Letter we report a detailed survey of the os- 
cillatory motion of rising bubbles in a wormlike micellar 
fluid. To study the role of the aggregation kinetics of 




FIG. 1: Cusp shape change of an oscillating bubble rising 
through a 15mM CPCl/NaSal (weight fraction tp = 0.5%) 
solution at T = 37.5° C. The scale at left is marked in cen- 
timeters. Interval between pictures: 0.05 s. 



the micelles, both concentration and fluid temperature 
were varied over a wide range, unlike previous studies 
[Til Ip3|. As these parameters change, different micellar 
architectures are possible - from short linear or branched 
micelles to crosslinked networks ^3, 13 01} ■ Four dif- 
ferent dynamics were seen: Newtonian behavior at high 
temperatures, standard polymeric behavior, and two dis- 
tinct oscillating responses occurring in different concen- 
tration ranges (Fig|2J). We performed steady rheology 
experiments to identify the fluid microstate, and found 
that transitions in the equilibrium structure match tran- 
sitions in bubble dynamics. Critical temperature bounds 
were also found, which can be interpreted as a minimum 
length of micelle required for oscillations to occur, so that 
for certain concentrations the fluid may be tuned with 
temperature to make bubbles either oscillate, rise with a 
stable cusp, or rise as bubbles in a Newtonian fluid. 

Non-transient oscillating bubbles were first observed in 
aqueous solutions of cetyltrimethylammonium bromide 
(CTAB) and sodium salicylate (NaSal) here we use 
another familiar system, cetylpyridinium chloride (CPC1) 
and NaSal, with the fixed ratio [NaSal]/[CPCl] = 1 (ex- 
cept for the NaCl experiments described below). Our 
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FIG. 2: Temperature and concentration phase diagram for 
the dynamics of rising bubbles in equimolar CPCl/NaSal, 
showing two distinct regions of oscillating behavior (shaded) 
labelled as I and II. The thick sloped line marking the tem- 
perature boundary for type I oscillations is an isoline of Eqn. 
1. 
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choice for this ratio gives a solution of flexible cylindrical 
micelles analogous to polymer chains [f| . All experimen- 
tal results are for a single rising bubble, with volumes 
ranging from 14 mm 3 to 110 mm 3 . The bubble is injected 
into the fluid at the bottom of a temperature controlled 
plastic cylinder (31 cm height, 5 cm diameter) using a 
syringe with a long stainless steel tube (inner diameter 
of either 1.0 mm or 1.5 mm). The ratio of the horizon- 
tal diameter of the bubble d to the cylinder diameter D 
is d/D < 0.14, and we have also checked that bubbles 
oscillate for d/D ~ 0.02. Typical sizes of the Reynolds 
number (inertia) are Re ~ 10~ 2 — 5, while the Debo- 
rah number (elasticity) ranged from De ~ 1 — 500, sim- 
ilar to values seen for oscillating bubbles and spheres in 
CTAB/NaSal [Hill- Rheological data were taken with 
a Rheometrics RFS-III controlled rate of strain rheome- 
ter with circulating fluid temperature bath, using a stain- 
less steel Couette geometry. Video images were made 
with a Kodak high speed CCD camera. 

We have studied concentrations from 4 to 40 mM 
(weight fractions 0.13% < <p < 1.3%). From 5 to 15 
mM, bubbles have a cusp which oscillates in length and 
changes shape (Fig.l). While rising, the cusp length- 
ens (frames 1-4 in Fig.l) during which the velocity as 
measured at the top of the bubble is nearly constant. 
At the apex of the extension, the tail abruptly retracts 
and the bubble jumps upward (frames 4-5). After this 
it slows to a nearly constant velocity until the cycle re- 
peats. Typically, bubbles begin to oscillate within 10 
seconds of their formation, with a similar time between 
jumps. These clearly visible oscillations, which we call 
type I, are apparently not a transient effect; we have ob- 
served their persistence for rise distances over one meter, 
during which there were more than 30 oscillations in ~ 
35 s (for 8 mM CPCl/NaSal), and velocities more than 
doubled during a jump. This is the typical type I behav- 
ior in the temperature ranges shown in FigEl consistent 



FIG. 3: Rheology of equimolar CPCl/NaSal at 30°C: a) ef- 
fective viscosity versus shear rate for 10, 20, 30, and 40 mM; 
b) zero shear viscosity, rjo as a function of concentration. The 
shaded regions mark the concentration ranges of bubble os- 
cillations, and the sloped line starting at 3 mM corresponds 
to r/ ~ V3 5 ' 8 . 



with the oscillations observed in CTAB/NaSal [Tlj . 

Above a critical temperature in the type I range 
(Fig|2J , bubbles have a sharp cusp which does not change 
in length ( "steady cusp" in Fig|2J and velocities smoothly 
reach a steady state; the oscillatory instability vanishes. 
By 60°C, the solutions appear Newtonian ("no cusp" 
in FigFJJ): large enough bubbles are ellipsoidal and un- 
dergo the well-known side-to-side oscillations In this 
temperature range, the micelles must be predominantly 
spherical or short rigid rods [2l| . 

Upon increasing the concentration from 15 mM to 16 
mM (ip — 0.5% to 0.53%), all oscillations cease. Bubbles 
rise steadily with a stable cusp for concentrations up to 
20 mM; there are no temperatures at which oscillations 
occur (Fig|5J). It is striking that a rising bubble would be 
so affected by so slight a concentration change. Note the 
transition temperature from polymeric to Newtonian be- 
havior reaches a maximum at this concentration (Fig[2J). 

More surprising is the re-emergence of oscillations for 
concentrations from 25 mM to 40 mM (0.8% < <p < 
1.3%). Yet here the oscillations, which we call type II, 
are visibly different from type I oscillations; the shape 
change involves the entire bubble, whereas in type I it 
seems restricted to the tail. A type II oscillation begins 
with a constriction in width near the top, while the whole 
bubble lengthens. This constriction then travels down- 
ward, as if the bubble were squeezing through a hoop. 
Defining w max to be the maximum (relaxed) width (at 
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its waist), w m in to be its most constricted width, and 
Aw = w max - w min , we found Aw/w max = 0.13 in 
type II fluids and 0.04 in type I. Length extensions, how- 
ever, were Al/l m i n = 0.25 in type II and 0.27 in type I; 
for bubbles rising in CTAB/NaSal, Al/l mm = 0.26 [nj. 
This previous study revealed only one type of oscillation, 
which from the shape dynamics we identify as type I |llj | . 

Although wormlike micelles are comprised of surfac- 
tants, it is unlikely that surface tension plays a role in the 
bubble oscillations, since falling rigid spheres also oscil- 
late |l2j |. Furthermore, our observation of two oscillation 
types suggests that the viscoelastic character of the bulk 
fluid is changing with concentration. We address this 
with rheology measurements, controlling the shear strain 
rate 7 and recording the effective viscosity, rj c ff — <r xy / A f 1 
where a xy is the steady shear stress. Transient tests were 
also performed to ensure that steady state was achieved. 
Shown in FigOJt is 7/off vs. 7 at 30°C, for fluids in the os- 
cillating and non-oscillating concentration ranges. Most 
fluids are shear thinning - r) c g decreases at high enough 
7 - however fluids below 5 mM shear thicken (r] e g in- 
creases above 770), typically requiring ~ 100 seconds to 
reach steady state [{j. Bubbles tested in fluids below 5 
mM showed no oscillations or steady cusps. While there 
is a noticeable difference in the rheology upon increas- 
ing from 20 to 30 mM, 30 and 40 mM (type II region) 
are strikingly similar; they also have a constant viscos- 
ity over a broader 7 range than the other fluids. More 
interestingly, the zero shear viscosity 770 (defined as the 
plateau value at low 7) decreases with concentration. 

A more complete study of the dependence of r/o on ip in- 
dicates transitions in fluid microstructure (Figl^). The 
transition at low concentration to rapid 770 growth marks 
the overlap concentration ip* ~ 3 mM, below which is the 
dilute regime. Above <p* is the semi-dilute state, in which 
micelles exist as individual entangled worms 22]. Here 
7/0 ~ ip 5 - 8 , close to the value 5.4 associated with stress 
relaxation by reptation 0, [2j| . The semi-dilute regime 
ends near tp ~ 12 — 15 mM, followed by an extraordinary 
decrease of 770 with concentration, continuing to tp ~ 30 
mM, after which 770 varies weakly. 

The dramatic change in the concentration dependence 
of 770 indicates a transition in the equilibrium fluid struc- 
ture, coinciding with the loss of type I oscillations. It may 
be that the entangled micelles in the semi-dilute range 
have begun to fuse at entanglement points, forming a 
crosslinked network |18l I23L l24l l25j. The junction nodes 
are free to slide along the micelles in such a formation, 
which would account for the decrease in viscosity |18| . If 
the ratio of crosslinks to entanglement points grows as 
770 decreases, then the new state would be fully formed 
where 770 stabilizes near ~ 30 mM (1%), the start of type 
II oscillations. Note that a crosslinked network state has 
also been proposed for CTAB/NaSal at <d ~ 1% |26|. 

Transitions in micellar morphology would naturally 
lead to transitions in mechanisms for stress relaxation, 
which should be observable in the stress field around the 
rising bubble. These transitions are made evident with 
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FIG. 4: Birefringent images of the wake behind a bubble 
rising in CPCl/NaSal at T = 24°C: a) 10 mM ; b) 20 mM ; 
c) 35 mM. Each image is 6.5 cm high. Reynolds and Deborah 
numbers for each image are: a) Re ~ 4.72, De ~ 250; b) 
Re ~ 0.02, De ~ 6; c) Re ~ 1.1, De ~ 1.8. 



birefringent visualization, shown in Fig^ for the three 
concentration regimes. Optical birefringence is a well- 
known technique for visualizing stress in non-Newtonian 
fluids [2]J , especially effective for wormlike micellar fluids 
|2£jj . FigQJ, (10 niM) shows the localization of stress in 
the wake of a type I oscillating bubble. This birefringent 
tail mirrors the dynamics of the bubble's cusp (Fig^l, 
and when it retracts, the birefringent tail suddenly dis- 
perses to the sides. The flow in the bubble wake is nearly 
in uniaxial extension |29|. and the birefringent tail sug- 
gests a strong relation to a thin filament whose rupture 
[l(j may be related to the mechanism for oscillation. At 
20 mM a shorter tail is seen (FigQJ)), whose length re- 
mains constant (like the bubble tail). Fig^J; corresponds 
to type II oscillations (35 mM), in which three equally 
spaced birefringent bands indicate a broader distribution 
of stress. This pattern is essentially unaltered during an 
oscillation. Note that De decreases from type I to type 
II, with an intermediate value for non-oscillating bub- 
bles (Fig. 4 caption), contradicting the implicit critical 
De condition for oscillations in |l2l ]. 

The length distribution of the individual worms in the 
type I region (semi-dilute) should depend on temperature 
and concentration. In equilibrium, the average length of 
a wormlike micelle Ln can be described with the mean 
field approximation |2ll3Ll2l|: 

L^T)~^>e E ' 2kT (1) 

where <fi is the total amphiphile concentration, k is Boltz- 
mann's constant, T is temperature, and E is the scis- 
sion energy required to break one wormlike micelle into 
two. Estimates based on light scattering measurements 
for E arc on the order of lOkT 0, 0, H(j - much lower 
than the covalent bonds of polymers, yet large enough 
for some micelles to reach appreciable lengths against 
thermal fluctuations 0. The flow near the bubble may 
increase the equilibrium length, and there is evidence 
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that shear-induced structures (SIS) much larger than in- 
dividual micelles form in such flow in both CTAB/NaSal 
and CPCl/NaSal We assume simply 

that for the type I region, Lq is the relevant quantity 
determining if the bubble will oscillate. Specifically, we 
attribute the oscillations to a breaking instability of the 
elongated micelles or SIS in the wake, occurring only if Lq 
exceeds some critical value. Using Eqn.l, the solid line in 
FigEJcorresponds to a scission energy of E = 1.01 x 10~ 19 
Joules, about 24kT. 

Type II oscillations have a much simpler dependence 
on temperature; they occur only for T < T c = 36° C 
(FigJ^J. If crosslinks have indeed formed for T < T c , 
then the mean field length is no longer an appropriate 
quantity. The transition temperature T c suggests a crit- 
ical energy condition (E c ~ 4.3 x 10~ 21 Joules), which 
may correspond to the energy difference between junc- 
tions and endcaps. Thus for T > T c , endcaps would 
dominate and the fluid is comprised of individual mi- 
celles which can entangle like polymers. This possible 
explanation is consistent with both the rheology and the 
rising bubble dynamics. 

The two different oscillations we have observed appear 
linked to two different microstructures. We extended our 
study to include the well characterized ternary system 
CPCl-NaSal diluted in concentrated NaCl, with a ratio 

[NaSal]/[CPCll = 0.5 [HHH^. We tested several con- 
centrations |36j of these fluids for jumping bubbles in 
our apparatus, at T = 30°C, a temperature central to 
both type I and II oscillations (Fig|2Jl. We observed no 
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